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ABSTRACT 

When  gases  are  admitted  rapidly  to  a  vessel  the  temperature 
may  rise  several  hundred  degrees  in  some  circumstances.  A  simple 
formula  involving  the  specific  heat  ratio,  thermal  conductivity 
and  molecular  weight  of  a  gas  is  derived  which  successfully 
describes  the  maximum  value  of  this  temperature  pulse  as  a 
function  of  temperature  and  pressure  for  a  variety  of  gases. 
Safety  applications  and  corrections  to  measurements  of  ignition 
limits  and  ignition  delay  times  are  discussed. 
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I .  INTRODUCTION 


When  gases  are  admitted  rapidly  into  an  evacuated  vessel  the 
temperature  at  the  center  of  the  vessel  can  rise  dramatically 
hundreds  of  degrees  above  the  temperature  of  its  surroundings. 

If  this  is  not  realized  and  isothermal  conditions  are  assumed 
to  hold,  large  experimental  errors  can  be  made  in  the  determination 
of  ignition  limits.  Figure  1  shows  the  explosion  limits  of  a 
hydrogen-oxygen  mixture  with  large  discrepancies  in  the  upper 
limits  as  determined  by  static  heating  (dashed  line)  and  rapid 
admission  of  gas  to  a  pre-heated  vessel  through  a  solenoid 

i  * 

valve  (solid  line) 

If  the  gases  form  a  combustible  mixture  and  no  provision 
for  explosion  has  been  made  the  consequences  could  be  disastrous. 
Industrial  accidents  have  occured  during  the  transfer  of  potentially 
explosive  gases  from  one  vessel  to  another.  Sometimes  such 
accidents  are  blamed  on  the  buildup  of  static  charge  in  spite  of 
the  fact  that  adequate  electrical  grounds  have  been  provided.  A 
more  probable  cause  is  the  phenomenon  described  in  this  report. 

In  summary,  both  the  combustion  chemist  and  the  safety 
engineer  can  benefit  by  a  better  understanding  of  what  is  happening. 
It  is  the  purpose  of  the  present  report  to  contribute  to  this 
understanding. 

II.  THE  PHENOMENON 

The  phenomenon  has  been  carefully  investigated  by  Fine,  Gray 
2-5 

and  MacKinven  in  the  course  of  their  work  on  spontaneous  ignition 

Figure  2  gives  a  sample  of  their  measurements  of  the  maximum  tempera 

ture  increase  At  above  the  vessel  walls  reached  at  the  center  of  an 

3  3  o 

evacuated  spherical  vessel  (1.222  x  10  cm  )  at  190  C  as  a  function 


*Referenoes  are  listed  on  page  28. 
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P (TORR) 


T  (°C) 

Figure  1.  Explosion  Limits  of  a  Hydrogen-oxgyen  Mixture  as  Determined 
bv  Static  Heating  (dashed  line)  and  Rapid  Admission  of  Gas 
to  a  Pre-Heated  Vessel  (Solid  Line) 
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of  the  final  pressure  in  the  vessel  for  various  gases.  Smaller 

temperature  rises  occur  closer  to  the  vessel  walls  with  a 

parabolic  profile  along  a  vessel  diameter.  Generally  speaking, 

the  greater  the  pressure  rise  the  greater  the  maximum  of  the 

transient  temperature  pulse  with  a  larger  effect  for  monatomic 

gases  than  for  diatonic  gases  and  with  small  changes  for  heavy 

polyatomic  gases.  Within  a  monatomic  or  diatonic  series  however 

the  lighter  atoms  usually  have  a  smaller  effect.  Figure  3  shows 

2 

the  same  relations  for  a  vessel  at  room  temperature  .  Relations 

within  a  series  are  brought  out  more  clearly  by  Figure  4  where 

&T  is  plotted  against  the  ratio  of  specific  heats,  for  the  same 

experimental  conditions  as  Figure  2  but  with  the  final  pressure 

2 

only  equal  to  10  torr  .  Again  it  is  clear  that  the  heavier  atoms 
show  a  larger  effect  in  the  monatomic  series,  however  in  the 
diatomic  series  oxygen  shows  a  smaller  effect  than  nitrogen  in 
spite  of  the  fact  that  it  is  more  massive.  Still  both  show  a 
larger  effect  than  hydrogen.  The  fact  that  there  are  other  factors 
at  work  besides  the  molecular  weight  and  y  of  a  molecule  is  made  clear 
by  the  observation  that  AT  for  argon  increases  with  an  increase  in 
vessel  temperature  while  AT  for  neon  and  helium  decrease.  Figure  5 
shows  the  same  relations  as  Figure  4  under  the  same  conditions  except 
that  the  vessel  is  at  room  temperature  instead  of  190°C. 

All  of  these  gases  were  studied  under  identical  conditions  and 

were  admitted  to  the  vessel  through  a  solenoid  valve  with  the 

temperature  pulses  for  a  final  pressure  of  10  torr  reaching  their 

maximum  values  in  about  0.1  or  0.2  seconds  and  decaying  almost  as 

rapidly.  The  temperature  changes  took  place  throughout  the  volume 

but  the  changes  at  the  center  were  always  greatest.  Cooling  effects 

were  also  observed  on  expansion  of  gases  from  a  vessel  at  higher  than 

2 

atmospheric  pressure  .  In  summary,  the  size  of  the  effect  is  a 
function  of  the  properties  of  the  gases  themselves,  the  temperature 
and  pressure  employed  and  the  speed  of  transfer.  In  the  next  section 
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P(TORR) 

Figure  2.  Maximum  Temperature  Increase  aT  at  the  Center  of  an  Evacuated 
Spherical  Vessel  (Wall  Temperature  190°C)  as  a  Function  of 
Final  Pressure  for  Various  Gases) 
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P(TORR) 

Figure  3.  Maximum  Temperature  Increases  AT  at  the  Center  of  an  Evacuated 
Spherical  Vessel  at  Room  Temperature  as  a  Function  of  Final 
Pressure  for  Various  Gases. 
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Figure  4.  Maximum  Temperature  Increase  AT  at  the  Center  of  an  Evacuated  Spherical  Vessel  (wall 

temperature  190°C)  for  a  Final  Pressure  of  10  Torr  as  a  Function  of  the  Ratio  of  Specific 
Heats,  for  Various  Gases:  1,  argon;  2,  neon;  3,  helium;  4,  nitrogen;  5,  oxygen;  6,  carbon 
monoxide;  7,  hydrogen;  8,  methane;  9,  ethane;  10,  propane;  11,  diethyl  ether. 
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Figure  5.  Maximum  Temperature  Increase  AT  at  the  Center  of  an  Evacuated  Spherical  Vessel  at  Room 
Temperature  for  a  Final  Pressure  of  10  Torr  as  a  Function  of  the  Ratio  of  Specific  Heats 
for  Various  Gases:  1,  argon;  2,  neon;  3,  helium;  4,  nitrogen;  5,  oxygen;  6,  carbon  monoxide 
7,  hydrogen;  8,  methane;  9,  ethane;  10,  propane;  11,  diethyl  ether. 
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of  this  report  we  shall  develop  a  simple  formula  to  describe  the 
observed  behavior. 


III.  THEORY  OF  THE  PHENOMENON 
A.  A  General  Formula 

When  a  vessel  of  gas  at  higher  pressure  is  connected  to  one 
at  lower  pressure  we  might  think  of  the  gas  in  the  first  vessel  as 
having  a  potential  energy  relative  to  the  final  state  in  which  both 
are  at  the  same  pressure.  As  soon  as  the  valve  connecting  the  two 
vessels  is  opened  the  gas  in  the  first  vessel  converts  its  potential 
energy  into  kinetic  energy  and  rushes  into  the  second  vessel  in  the 
form  of  a  jet.  The  kinetic  energy  of  this  jet  is  then  converted  into 
internal  energy  of  the  gas  which  has  been  admitted  by  collisions  with 
the  vessel  walls  and  with  other  molecules.  This  manifests  itself  as 
a  temperature  rise  which  reaches  its  maximum  value  AT  at  the  center 
of  the  vessel  where  heat  exchange  with  the  surroundings  is  slowest. 

Let  us  write  this  increase  in  internal  energy  as  (1/2)  fNkAT  where 

f  is  the  number  of  degrees  of  freedom  of  the  molecule,  N  is  the 

number  of  molecules,  k  is  Boltzmann's  constant  and  AT  is  the  maximum 

size  of  the  thermal  pulse  at  the  center.  The  energy  which  has  been 

admitted  into  the  vessel  is  in  the  form  of  n  moles  of  gas  with  kinetic 
2 

energy  (l/2)NmvB  ,  some  of  which  is  lost  to  the  walls  by  heat  transfer 
with  coefficient  K(erg  per  mole-deg) .  We  may  then  write 

(1/2) fNk  AT  -  (l/2)Nmvf  -  nK  AT  (1) 

Eu 

as  the  average  energy  balance.  Here  we  are  using  the  most  probable 

1/2 

speed  of  the  molecules  v  ■  (2RT  /M)  instead  of  the  root  mean 

1/2™  " 

square  speed  (3RTw/M)  in  terms  of  the  molecular  weight  M  and  the 
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temperature  of  the  vessel  wall  T  with  AT  ■  T  -  T  where  T  is  the 

w  w 

maximum  temperature  reached  at  the  center  of  the  vessel.  All 
except  the  heavy  gases  are  heated  to  the  temperature  of  the  vessel 
walls  in  the  entry  arm.  However,  propane  and  diethyl  ether  show 
slight  cooling  before  pulse  heating  since  they  are  not  fully  warmed 
on  entry  because  of  their  high  heat  capacities.  This  results  in  a 
lower  AT  than  expected.  Since  (1/2) f  ■  1/(Y-1),  Nm  ■  nM  and  NK  »  nR, 
we  obtain 

AT  -  (Y-l)Tw/(l  ♦  (K/R)(y-1)).  (2) 

Alternatively  we  may  consider  the  gas  admitted  first  to  be 
compressed  by  the  gas  admitted  later  and  use  the  first  law  of 
thermodynamics  for  an  ideal  gas  in  the  form 

Aq  -  cyAT  ♦  pAv  ■  cyAT  ♦  RAT  -  RTAp/p  (3) 

where  v  is  the  specific  volume  V/n  and  T  is  the  temperature  after  a 
given  amount  of  work  has  been  done.  If  we  take  Aq  »  -KaT, 
c  ♦  R  ■  c  :  ■  yR/(y-1)  and  Ap  ®  p  since  the  initial  pressure  in  the 
vessel  is  practically  zero,  we  find 

AT  -  T  -  Tw  -  (y-1)T/(y*  (K/R)(y-1))  (4) 

or  in  terms  of  T 

w 

AT  -  (Y-1)Tw/(1  ♦  (K/R)(y-1)  )  (5) 

as  in  equation  (2)  above.  If  the  vessel  already  contains  some  gas, 

Ap/p  would  be  an  appropriate  fraction  r  instead  of  unity.  Equation  (5) 
would  then  be 

AT  =  r(Y-l)Tw/(y  (K/R  -  r)(Y-l)  ).  (6) 
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This  can  also  be  derived  from  the  first  aodel  if  we  multiply 
2 

(l/2)Nmvffl  by  the  fractional  contribution  r  to  the  total  energy 
in  the  vessel.  In  other  words,  the  added  energy  must  be  shared 
between  the  molecules  added  to  the  vessel  and  those  already  inside. 

The  effect  of  r  less  than  unity  is  to  lower  AT. 

Both  derivations  are  crude  energy  balances  which  take  no 
account  of  the  complicated  changes  which  are  taking  place  in  time 
and  space.  However  our  only  aim  is  to  predict  the  maximum  temperature 
rise  at  the  center  of  the  vessel.  In  addition  we  will  treat  the  timing 
factor  in  an  empirical  fashion  leading  simply  to  useful  results. 

B.  The  Heat- loss  Factor 

Heat  losses  occur  to  the  vessel  walls  by  a  combination  of 

2 

convection  and  conduction.  Studies  have  been  made  which  suggest 

that  conduction  is  dominant  at  pressures  up  to  100  torr  at  least 

during  the  cooling  portion  of  the  temperature  pulse.  Let  us  write 

our  thermal  loss  coefficient  K  ■  XLt/n  where  *  is  the  coefficient  of 

thermal  conductivity  (erg/cm-sec-deg) ,  L  is  the  effective  loss  distance 

(cm)  and  r(sec)  is  the  time  for  the  gas  to  enter  and  the  temperature 

pulse  to  reach  its  maximum  at  the  center.  If  the  gas  is  admitted  very 

slowly,  K  becomes  very  large  and  AT  approaches  zero,  a  practical  way 

to  avoid  unwanted  temperature  pulses,  but  a  poor  way  to  determine 

delay  times  or  ignition  limits  in  a  reaction.  Since  n  ■  pV/(RTw) 

we  ha-'e  K/R  -  (L/VJx-rT^/p,  showing  that  losses  are  smaller  at  higher 

pressures  giving  larger  AT.  If  we  use  the  experimental  value  t*  .21 

3  3 

seconds  for  nitrogen  at  10  torr  final  pressure  in  a  1.222  x  10  cm 
spherical  vessel  at  room  temperature  which  produces  a  pulse  AT  ■  57°C, 
we  find  that  K/R  must  be  2.86  which  gives  a  value  of  L  ■  280  cm,  much 
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larger  than  the  vessel  diameter  of  about  10  cm.  In  the  context  of 

the  semi-empirical  approach  we  have  adopted  we  might  attribute  this 

to  the  effect  of  convection  on  the  heating  portion  of  the  temperature 

pulse.  Forced  convection  must  make  a  considerable  contribution  when 

a  jet  of  gas  rushes  into  an  evacuated  vessel  although  density 

differences  and  so  free  convection  may  be  unimportant.  In  the 

simple  model  we  have  developed  we  will  represent  convection  through 

the  factor  L  which  we  will  fix  at  the  value  above  in  all  our 

calculations.  Where  experimental  values  of  t  are  not  available 

we  will  calculate  them  from  the  appropriate  values  for  nitrogen 

1/2  2 

by  the  formula  T  =  TN  (M/28)  since  it  has  been  observed  that 
T  is  approximately  proportional  to  the  square  root  of  the  molecular 
weight.  It  is  likely  that  t  is  inversely  proportional  to  the  molecular 
speed  which  would  make  it  proportional  to  (M/Tw)  '  .  This  agrees 
approximately  with  the  observed  values  (Cf. Table  I)  and  enables  us 
to  calculate  T  as  a  function  of  temperature  also. 

IV.  COMPARISON  WITH  EXPERIMENT 

Table  I  gives  the  required  parameters  for  calculating  AT  for  a 
final  pressure  of  10  torr  for  a  variety  of  gases  and  compares  these 
calculated  values  with  experimental  values  taken  from  Figures  4  and  5. 
As  we  see,  if  we  omit  propane  and  diethyl  ether  which  show  an  initial 
cooling  before  heating  the  average  deviation  is  about  1.33°C  with  a 
maximum  deviation  of  7°C.  If  we  had  used  a  calculated  value  for  Helium 
at  300°K  instead  of  the  quoted  experimental  value,  the  7°C  deviation 
would  have  been  much  reduced.  This  is  remarkably  close  agreement 
considering  the  crudeness  of  the  model  and  the  accuracy  of  the 
experiment  and  encourages  us  to  compare  the  predictions  of  the  model 
with  experiments  as  a  function  of  pressure  (Figures  2  and  3) . 
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Table  I.  Required  Parameters  M  (gm/mole) ,A(erg/cm-sec-deg) , t  (sec)Y  ,  and  Experimental  and  Calculated 
Values  of  (°K)  for  Admission  of  Gas  to  a  Final  Pressure  of  10  Torr  at  Two  Temperatures  of 


the  Initially  Evacuated  Spherical 

M  Y300  Y463  X300Xl° 

Reaction  Vessel 

*465  10'4 

e 

T300 

T463 

AT300 

AT463 

A 

40 

1.67 

1.67 

0.18 

0.25 

.225 

.185 

77(82) 

87 

(88) 

Ne 

20 

1.67 

1.67 

0.49 

0.65 

.147 

.130 

56(57) 

48 

(55) 

He 

4 

1.67 

1.67 

1.50 

2.01 

.085 

.075 

43(36) 

32 

(33) 

N2 

28 

1.40 

1.40 

0.26 

0.37 

.210 

.165 

57(57) 

64 

(62) 

°2 

32 

1.40 

1.37 

0.27 

0.39 

(.220) 

(.177) 

50(53) 

58 

(57) 

CO 

28 

1.40 

1.40 

0.25 

0.36 

(.210) 

(.165) 

57(57) 

63 

(63) 

»2 

2 

1.40 

1.40 

1.82 

2.46 

(.056) 

(.044) 

40(41) 

35 

(40) 

ch4 

16 

1.32 

1.23 

0.34 

0.60 

(.159) 

(.125) 

45(50) 

45 

(44) 

C2H6 

30 

1.18 

1.12 

0.22 

0.45 

(.217) 

(.171) 

35(37) 

27 

(32) 

C3H8 

44 

1.16 

1.09 

0.18 

0.37 

(.260) 

(.206) 

21(35) 

16 

(27) 

Cc2h5)2° 

74 

1.09 

1.06 

0.15* 

0.32 

(.335) 

(.270) 

14(22) 

7 

(20) 

‘Extrapolated  ignoring  boiling  point  =  308°K.  Values  in  parentheses  were  calculated 


Table  11.  Experimental  and  Calculated  Values  of  AT(°K)  as  a 
Function  of  Pressure  for  an  Initially  Evacuated 
Reaction  Vessel  at  Two  Temperatures. 

T  =  300°K 

w 


pressure 

(torr) 

10 

20 

30  40 

50 

60 

80 

100 

argon 

77 

(82) 

no 

(117) 

130  150 

(135) (147) 

155 

(156) 

160 

(163) 

175 

(171) 

(175) 

neon 

56 

(57) 

80 

(88) 

105  125 

(109) (122) 

135 

(132) 

145 

(141) 

155 

(152) 

(160) 

helium 

43 

(36) 

65 

(62) 

80  105 

(80) (  94) 

115 

(105) 

125 

(114) 

(128) 

(138) 

nitrogen 

57 

(57) 

75 

(76) 

80  85 

(87)  (93) 

T  «  463°K 
w 

95 

(98) 

100 

(101) 

105 

(105) 

(108) 

pressure 

(torr) 

10 

20 

30  40 

50 

60 

80 

100 

argon 

87 

(88) 

140 

(137) 

165  180 

(168)  (190) 

200 

(207) 

210 

(218) 

230 

(237) 

245 

(248) 

neon 

48 

(55) 

100 
(  94) 

125  150 

(123) (144) 

165 

(161) 

180 

(176) 

200 

(196) 

210 

(212) 

helium 

32 

(33) 

70 

(61) 

100  125 

(  83) (102) 

130 

(118) 

150 

(131) 

170 

(150) 

185 

(170) 

nitrogen 

64 

(62) 

85 

(93) 

100  no 

(111) (123) 

120 

(132) 

130 

(139) 

135 

(148) 

(154) 

Calculated  values  are  in  parentheses.  Experimental  values  have  been 
estimated  from  Figures  2  through  5. 
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Table  II  compares  calculated  and  experimental  values  of  AT  as  a 
function  of  pressure  for  two  temperatures.  The  agreement  is  generally 
good  considering  the  model  and  the  experiment  and  might  have  been  better 
if  experimental  values  of  t  were  available  as  a  function  of  pressure. 
However,  the  values  of  t  measured  or  calculated  for  a  final  pressure  of 
10  torr  were  used  throughout  (as  was  a  constant  L  value)  and  the  K/R 
values  of  10  torr  were  taken  simply  as  inversely  proportional  to  the 
pressure,  probably  accounting  for  much  of  the  variance  observed. 

V.  DISCUSSION 

The  model  is  too  crude  to  warrant  many  refinements,  yet  it 

successfully  predicts  trends  including  the  observed  crossover  of  the 

nitrogen  curve  with  the  neon  and  helium  curves  as  a  function  of 

pressure  (Figures  2  and  3).  It  also  predicts  an  upper  limit  to  the 

temperature  increase  of  AT  =  (y-l)T  when  K  =  0.  When  T,  =  300°K, 
v  max  w  w 

AT  =  200°K  for  monatomic  gases,  12CPK  for  diatomic  gases,  and  so 

IHoX 

forth.  When  Tw  =  463°K  (190°c) ,  ATmax=  310°K  for  monatomic  gases  and 
185°k  for  diatomic  gases,  etc.  These  predictions  are  of  special 
interest  from  a  safety  viewpoint  and  slow  gas  admission  can  be  used  to 
eliminate  danger.  When  the  measurement  of  ignition  delays  and  ignition 
limits  is  of  interest  we  can  do  several  things  to  reduce  error.  If  a 
heavy  fuel  plus  oxygen  system  is  being  studied  the  oxygen  might  be 
admitted  to  the  reaction  vessel  before  it  is  heated  and  the  fuel 
admitted  rapidly  after  the  desired  temperature  has  been  reached  since  a 
heavy  molecule  (unlike  oxygen)  will  only  give  a  small  temperature  pulse. 
Alternatively,  if  we  desire  to  admit  the  oxygen  rapidly  into  a  vessel 
containing  the  fuel  we  can  get  an  estimate  of  the  temperature  pulse  by 
using  an  inert  gas  such  as  nitrogen  (or  an  inert  model  of  the  fuel  if 
the  fuel  were  to  be  transferred) .  Or  in  some  circumstances  we  might 
adopt  the  idea  of  two  vessels  connected  by  a  wide-bore  stopcock  with 
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diffusion  mixing  after  pressure  equilibration.  Finally,  if  premixed 
gases  are  admitted  to  a  vessel  effective  average  values  of  M,  y,  and  X 

can  be  used  to  calculate  AT  and  make  corrections. 

VI .  CONCLUSION 

We  have  derived  a  simple  formula  for  calculating  the  temperature 
increase  which  results  when  a  gas  is  admitted  rapidly  to  a  vessel.  In 
spite  of  its  simplicity  the  predictions  of  this  formula  agree  well  with 
experimentally  measured  values  as  a  function  of  temperature  and 
pressure  for  a  variety  of  gases.  This  enables  us  to  understand  the 
phenomenon  more  clearly  and  to  take  precautions  or  make  corrections 

as  required. 
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